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FOR Z ~ 2.2 LYMAN ALPHA EMITTERS f * tt 

Takuya Hashimoto 1 , Anne Verhamme 2 , Masami Ouchi 3 ’ 4 , Kazuhiro Shimasaku 1,6 , Daniel Schaerer 2 ’ 6 , 
Kimihiko Nakajima 2 , Takatoshi Shibuya 3 , Michael Rauch 7 , Yoshiaki Ono 3 ’ 4 , and Ryosuke Goto 1 . 

ABSTRACT 

We present the results of a Lya profile analysis of 12 Lya emitters (LAEs) at z ~ 2.2 with high- 
resolution Lya spectra. We find that all 12 objects have a Lya profile with the main peak redward of 
the systemic redshift defined by nebular lines, and five have a weak, secondary peak blueward of the 
systemic redshift (blue bump). The average velocity offset of the red main peak (the blue bump, if 
any) with respect to the systemic redshift is AuL ya ,r = 174±19 km s _1 (AriLya.b = — 316±45 km s _1 ), 
which is smaller than (comparable to) that of Lyman-break galaxies (LBGs). The outflow velocities 
inferred from metal absorption lines in three individual and one stacked spectra are com parable to 
those of LBGs. The uniform expanding shell model constructed bv iVerhamme et all (|2006ft reproduces 
not only the Lya profiles but also other observed quantities including the outflow velocity and the 
FWHM of nebular lines for the non-blue bump objects. On the other hand, the model predicts too 
high FWHMs of nebular lines for the blue bump objects, although this discrepancy may disappear 
if we introduce additional Lya photons produced by gravitational cooling. We show that the small 
AiiLya.r values of our sample can be explained by low neutral-hydrogen column densities of log(IVHi) 

= 18.9 cm -2 on average. This value is more than one order of magnitude lower than those of LBGs 
but is consistent with recent findings that LAEs have high ionization parameters and low Hi gas 
masses. This result suggests that low TVhi values, giving reduced numbers of resonant scattering of 
Lya photons, are the key to the strong Lya emission of LAEs. 

Subject headings: galaxies: high-redshift — galaxies: ISM - line: profiles — radiative transfer 


1. INTRODUCTION 

Lya emitters (LAEs) are objects commonly seen 
in both the local and high -2 universes with large 
Lya equivalent w idths, EW(Lyq)> 2 0 — 3 0 A(local: 
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Deharveng et all 120081 iGowie et all 1201 ll h igh- z : 

ad 1996 : iRhoads fe Malhotral 120011 
2010 1. Previous studies based on 


Spectral Energy Distributions (SEDs) have revealed 
that typical LAEs are young, low-mass ga laxies 
with a small dust content (e.g.. [Ni lsson et ahl 120071 
iGawiser et al.ll2007l iGuaita et al.1 120111 iNakaiima et al.1 
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I2012t iKusakabe et al.l 12015ft , although there are some 
evolve d LAEs with a mode rate mass and dust (IQno et al.l 
l2010bl lHagen et al.1 12014ft . Morphological studies of 
their UV continuum have shown that the galactic 
counterparts o f LAEs to be typically compact (e.g., 
iBond et all 12009ft and their typical size does not evolve 
with redshift dMalhotra et al.ll2012li . Furthermore, clus¬ 
tering analyses have revealed that LAEs ha ve the lowest 
dark matter halo masses at every redshift (|Ouchi et al.l 
120101 : IGuaita et ahll2010ft . These properties suggest that 
LAE is an important galaxy population as the b uil ding 
block candidates in the A CDM model (I Rauch et al.1 

12008ft . 

Given their importance in galaxy evolution, the Lya 
escape mechanism in LAEs is still poorly understood. 
Resonant scattering strongly extends the path-length 
of Lya photons through galactic gas and renders them 
prone to absorption by dust grains. On one hand, some 
observational studies at the local universe have proposed 
that outflows fac il itate the esca pe of Lya photons from 
galaxies (e.g JKunth et al.lll998ft as they re duce the num¬ 
ber of scattering . Lik ewise, others (e.g., IKornei et al.1 
120101 lAtek et ahl 12014ft have shown that the dust con¬ 
tent correlates with Lya emissivity. While these effects 
would certainly be at work, there has been no deci¬ 
sive conclusion (cf., ICassata et al.l [2015ft . On the other- 
hand, theoretical studies have computed the Lya radia¬ 
tion transfer (RT) through idealized spherically symmet¬ 
ric shells of homogeneous and isothermal neutral hydro¬ 
gen gas, especially in a form of an expanding shel l (e.g. , 
Zheng fc Miralda-Escudel l2002t IVerhamme et al.l 120061 
Diikstra fe Loebll2009tlKollmeier et al.ll2010ft . They have 
investigated how properties of the interstellar medium 
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(ISM) affect the Lya escape and emergent Lya profiles. 
The result is that the Lya RT is a complicated process al¬ 
tered by galactic outflows/inflows, the neutral hydrogen 
column density and dust content of the ISM, and the in¬ 
clination of the galaxy disk (e.g.. IV erhamm e et alJl2012t 
IZheng fe Wallacell20l4 IBehrens fe Braunll2014T) . One of 
the goals in these theoretical studies is to aid understand¬ 
ing the galaxy properties from observed Lya lines, and 
to identify the key factor for the Lya escape. 

To study the Lya RT and escape through close com¬ 
parisons of observed and modeled Lya lines, it is im¬ 
portant to obtain spectral lines other than the Lya line. 
The central wavelength and the width of nebular lines 
(e.g., Ha and [Om]) tell us the galaxy’s systemic red- 
shift and internal velocity. The blue-shift of interstellar 
(IS) absorption lines with respect to the systemic red- 
shift gives the galactic outflow velocity, and the width 
of the IS lines can be interpreted as the sum of thermal 
and macroscopic (rotation and turbulence) velocities of 
the outflowing gas. These lines can help us to disentangle 
the complicated Lya RT and understand the Lya escape. 

However, due to the typical faintness of LAEs, it is 
only recently that these additional lines have been suc¬ 
cessfully detect ed in narrow-band selected LAEs (nebu¬ 
lar li n es: e.g.. IM cLinde n et all 120111 : iFinkelstein et al.1 
201 It lH ashimo to et al. 20131. IS absorp tion lines: 


Hashimoto et al.l 120131 : IShibuva et all I2014bh . Thus, in 


contrast to LBGs whose Lya profiles have been closely 
compared with Lva RT models fe.g.. IVe rhamm e et all 
120081 iKulas et al.l l2012t iChristensen et al.l 120121 1. there 
are only a few studies that h ave performed Lya pro¬ 
file comparisons of LAEs fe.g.. iChonis et ahll2013f) . Re¬ 
cent simultaneous detections of Lya and nebular emis¬ 
sion lines have statistically confirmed that the Lya pro¬ 
files of LAEs are asymmetric with a red main peak red- 
shifte d with respect to the systemic Aiii,y n r , > 0 km s -1 
(e.g., IShibuva et al.ll2014bl : iSong et alTl2014l : lErb et all 
1201411 . Likewise, IS absorption studies in LAEs have 
shown that they are blue-shifted with r espect to the sys¬ 
temic by |A^ a hsl ~ 100— 200 km s -1 ([Hashimoto et al.l 
12013 : Shibuvaet al.||201 4blf . which i s comparable to those 
of LBGs fe.g.. IP ettini et al.l 1200 11: Shaolev et al.l [20031 : 


iSteidel et al.ll201CtlKulas et al.ll2012 1. These results sug¬ 
gest that LAEs do have outflows and motivate us to apply 
expanding shell models to LAEs. 

To examine Lya escape mechanisms in LAEs through 
detailed Lya modeling, we focus on the small A-CLya.r 
of LAEs, AcLya.r — 200 km s -1 , compared to those of 
LBGs, AuT,vf>.r ^ 400 km s -1 (e.g., ISteidel et all 120101 : 
IKulas et al.ll2012li . with similar physical quantities such 
as stellar mass, star formation rate (SFR), and velocity 
dispersion (lHashimoto et~aJll2013l: Shibuva et al.l[2014bl: 


ISong et all 2014: Erb et all 2014 L lHashimoto et al l 

(|2013f ) and IShibuva et al.l ( 2014H 1 have also shown that 
LAEs have comparable outflow velocities, measured from 
IS absorption lines, to those of LBGs. These re¬ 
sults imply that a definitive difference between LAEs 
and LBGs in velocity proper ties is AtiL ya ,r. In addi¬ 
tion, lHashimoto et al.l ( 2013 1 have demonstrated that 
EW(Lya) anti-correlates with At>r, yn. r using a large sam¬ 
ple of LAEs an d LBGs (see also IShibuva et al.l l2014bl : 
lErb et all I2014H . Therefore, understanding the reason 
why LAEs have small AwLya,r, through detailed Lya 


modeling should shed light on the Lya RT and Lya es¬ 
cape mechanisms in LAEs. 

According to the theoretical studies, there are several 
possible expla nations for a smal l At>r, v r> r : a high outflow 
velocity (e.g., IVerhamme et al.l 12006 :1. a very low neu¬ 
tral hydrogen column densit y (IVhi) of the ISM (e.g., 
IVerhamme et al.l 120061 I2015H . an inhomogeneous ISM 
with a covering fraction (CF) below unity, where CF 
is defined as the fraction of sightlines which are opti¬ 
cally thick to the Lya radiation, i.e., gas with holes 
fe.g.. IBehrens et all 12014 IVerhamme et al.l 120151 1. and 
a clumpy ISM with a low covering factor, / c , which 
is defined as the average number of clouds intersected 
by a random line of si ght (e.g., lHansen fc Ohl 120061 : 
iDiikstra fc Kramerll2012l : TLaursen et al.ll2013fi 

In this work, we focus on applying the uniform ex- 
pa nding shell model based on a 3D Lya RT construc ted 
bv IVerhamme et ahl (12006T ) and lSchaerer et all (120111 ). to 
12 LAEs whose Lya and nebular emission lines (e.g., 
Ha, O m) have been obs e rved at a high spectral res¬ 
olution (lHashimoto et al.1 120131 : INakaiima et al.l 120131 : 
IShibuva et al.l l2014blf . With the systemic redshifts and 
the full width half maximums (FWHM) determined from 
nebular emission lines, the stellar dust extinction derived 
from SED fitting, and the galactic outflow velocities in¬ 
ferred from LIS absorption lines, we first statistically ex¬ 
amine how well the model can r eproduce the Lya pro¬ 
files and other observables fcf.. IVerhamme et al.l 120081 : 
IKulas et al .11 20 id iGhonis et al.ll2013h . After demonstrat- 

ing the validity of the model, we securely derive physical 
quantities such as Aqii and discuss the origin of the small 
Ar/Lya,r and implications for the Lya escape in LAEs. 
Possible other scenarios mentioned above for the small 
Ai/Lya,r are also qualitatively discussed. 

This paper is organized as follows. We describe our 
spectroscopy observations in Section [2J and discuss pro¬ 
files of Lya and nebular emission lines in Section [3] We 
apply the uniform expanding shell model to our data and 
show comparisons with observables in Section [4j Discus¬ 
sion on the blue bumps as well as the origin of the small 
Lya velocity offsets are given in Section [5l followed by 
conclusions in Section [Gl 

Through out this paper, mag nitudes are given in the 
AB system (|Oke fc Gunnlll983lf , and we assume a ACDM 
cosmology with H m = 0.3, Ha = 0.7 and H 0 = 70 km 
s -1 Mpc -1 . 


2. DATA AND OBSERVATIONS 

Our initial sample of objects are taken from large 
z ~ 2.2 LAE samples in the COSMOS field, the Chan¬ 
dra Deep Field South (CDFS), and the Subaru /XMM- 
Newt on Deep Survey (SXDS) (INakaiima et al.l 120121 . 
120131 : Nakajima et al. in prep.). These LAE samples are 
all based on Subaru/Suprime-Cam imaging observations 
with our custom made narrow band filter, NB387 (A c = 
3870A and FWHM = 94A). The LAEs have been selected 
by color criteria of B — NB387 and u* — NB387, satis¬ 
fying the condition that the rest frame photometric Lya 
EW (EW(Lya)photo) be larger than 30A. From these, 
we only use 12 LAEs whose Lya and nebular emission 
lines (e.g., Ha and [Om]) are both spectroscopically con¬ 
firmed. Among t h e 12 o bject s, 11 have been present ed in 
lHashimoto et al.l (120131) and IShibuva et al.l (l2014bD . We 
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add one new LAE with EW(Lya) p hoto ~ 280A whose 
detailed properties will be discussed in Hashimoto et al. 
in prep. 

In this section, we briefly summarize our near-infrared 
spectroscopy ( 92.11) . optical spectroscopy ( 92.2[) , and the 
contamination of AGNs in the sample ( 112.311 . 

2.1. Near-Infrared, Specctroscopy 

In order to detect nebular emission lines, we performed 
three near-infrared observations with Magcllan/MMIRS 
(PI: M. Ouchi), Keck/NIRSPEC (PI: K. Nakajima), and 
Subaru/FMOS (PI: K. Nakajima). Canonical spectral 
resolutions for our observation settings are R ~ 1120, 
1500, and 2200 for MMIRS, NIRSPEC, and FMOS, re¬ 
spectively. 

Details of the observation and data reduction proce¬ 
dures for MMIRS and N IRS PEC have been p r esent ed in 
lHashimoto et all (|2013l) and iNakaiima et all (|2013l h re¬ 
spectively. Briefly, two CDFS objects, CDFS-3865 and 
CDFS-6482, were observed with MMIRS using the HK 
grism covering 2.254 — 2.45 /im, resulting in success¬ 
ful Ha and [Oin]AA 4959 5007 detections. A follow-up 
observation was carried out for CDFS-3865 with NIR¬ 
SPEC. The [On] A3727 line was additionally detected 
with the J band (1.15 — 1.36 fin i). Four COSMOS 
objects, COSMOS-08501, COSMOS-13636, COSMOS- 
30679, and COSMOS-43982, were observed with NIR¬ 
SPEC and its K band (2.2 — 2.43 /rm), resulting in Ha 
line detections. The [Om] A5007 line was also detected 
from COSMOS-30679 using the H band (1.48 — 1.76 /jm). 

The data from FMOS will be presented in Nakajima et 
al. (2015, in prep). Its spectral coverage is 0.9 — 1.8 fim. 
We detected [Om] line(s) in eight objects: COSMOS- 
08357, COSMOS-12805, COSMOS-13138, COSMOS- 
13636, COSMOS-38380, COSMOS-43982, SXDS-10600, 
and SXDS-10942. 

2.2. Optical Spectroscopy 

In order to detect Lya and metal absorption lines, 
we carried out several observations with Magellan/MagE 
(PI: M. Rauch) and Keck/LRIS (PI: M. Ouchi). The 
spectral resolutions for our observations were R ~ 4100 
and ~ 1100 for MagE and LRIS, respectively. The slit 
was positioned on the Lya centroids in the NB387 im¬ 
ages. 

Details of the observation and data reduction pro¬ 
cedures for MagE and LRI S have been p resente d in 
lHashimoto et al.l (1201311 and iShibuva et al.l (l2014bfl , re¬ 
spectively, except for COS MOS-0 8501. First, we describe 
this new object in detail ( 92.2.11) and then giver a brief 
summary for the rest of the sample ( 92.2.211 . 

2.2.1. Optical Spectroscopy for COSMOS-08501 

The MagE observations were carried out for COSMOS- 
08501 on 2012 February and 2013 December. We ob¬ 
tained 3 x 3000 s and 1 x 3000 s exposure times during 
each run, resulting in a 12000 s total integration time. 
Spectroscopic standard stars, dome flats, and Xenon 
flash lamp flats, were obtained on each night for calibra¬ 
tions. On these nights, the typical seeing sizes were 1".0. 
The slit width was l."0 for both runs, corresponding to 
R ~ 4100. The spectra were reduced with IDL based 
pipeline, MagEJtEDUCE, constructed by G. Becker (see 


also iKelsonl l2003h . This pipeline processes raw frames, 
performing wavelength calibration and optimal sky sub¬ 
traction, and extracts ID spectra. Each of these reduced 
frames was then combined to form our final calibrated 
spectrum. From this, the Lya line was identified above 
the 3 a noise of the local continuum. 

2.2.2. Optical Spectroscopy for the Rest of the Sample 

CDFS-3865, CDFS-6482, and COSMOS-30679 were 
observed with MagE; COSMOS-08357, COSMOS-12085, 
COSMOS-13138, COSMOS-38380, SXDS-10600, and 
SXDS-10942 were observed with LRIS; and finally, 
COSMOS-13636 and COSMOS-43982 were observed 
with both spectrographs. We identified the Lya 
line in all objects. In addition, we detected sev¬ 
eral metal absorption lines (e.g., Si II A1260 and Civ 
A1548 lines) in a stacked MagE spectrum of CDFS- 
3865, CD FS-6482, COSMOS -13636, and COSMOS- 
30679 (lHashimoto et al.l 1201311 as well as in individual 
LRIS spectra of COSMOS-12805, C OSMOS-13636, and 
SXDS-10600 (IShibuva et al.H2014bll . 

A summary of our observations is listed in Tabic [TJ and 
our Lya and nebular emission line profiles are shown in 
Figure 1. 

2.3. AGNs in the Sample 

The prese nce of AGNs in the Ma gE objects has been 
exami ned in lHashimoto et all (|2013f) and lNakaiima et~ah 
( 2013I1. and those of the LRIS objects in IShibuva et al. 

In short, for the MagE objects, we inspected it in three 
ways. We first compared the sky coordinates of the ob¬ 
jects with those in very deep archival X-ray and radio 
catalogs. Then we checked for the presence of high ion¬ 
ization state lines such as Civ A 1549 and He ii A1640 
lines in the spe ctra. Fin ally , we app lied the BPT diag¬ 
nostic diagram dBaldwin et al.lIT98lh to the objects. No 
AGN activity is seen except for COSMOS-43982 whose 
high [Nil] /Ha line ratio is consistent with that of an 
AGN. 

On the other hand, due to the lack of Ha or [Nil] A6568 
data, we were only able to use the two forms of investi¬ 
gation for the LRIS objects. Of these, only COSMOS- 
43982 showed clear detection of the Civ A 1549 line in 
its optical spectrum. 

In summary, we have ruled out AGN activity in all but 
COSMOS-43982. 

3. OBSERVATIONAL RESULTS 

3.1. Line Center and FWHM Measurements for 
Nebular Emission Lines 

Line center (i.e., redshift) and FWHM measurements 
of nebular emission lines are crucial for a detailed model¬ 
ing of the Lya line, since they encode information on the 
intrinsic (i.e., before being affected by radiative trans¬ 
fer) Lya redshift and FWHM. In order to obtain these 
parameters and their uncertainties, we apply a Monte 
Carlo technique as follows. First, for each line of each 
object, we measure the ler noise of the local continuum. 
Then we create 10 * * 3 fake spectra by perturbing the flux 
at each wavelength of the true spectrum by the measured 
la error dKulas et al.ll2012tlChonis et al. 1120 131) . For each 
fake spectrum, the wavelength at the highest flux peak 
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Fig. 1.— The upper two rows of panels show the Lya lines (black) obtained with MagE and the corresponding Ha lines (red), while 
the lower two rows of panels show the Lya lines obtained with LRIS (black) and the [OIII] lines (red). Blue arrows denote blue bumps. 
All spectra are scaled to the wavelength range from —1500 to +1500 km s —1 . Yellow and green segments indicate the peak flux positions 
derived from a symmetric Gaussian and a Monte Carlo technique, respectively. 
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TABLE 1 

Summary of the Observations 


Object 

a(J2000) 

<5( J2000) 

EW (Lya)photo 

L(Lya) 

NIR obs. 

opt. obs. 

Source* 2 




(A) 

(10 42 erg s _1 ) 




(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

CDFS-3865 

03:32:32.31 

-28:00:52.20 

64 + 29 

29.8 + 4.9 

NIRSPEC (J) 
MMIRS (HK) 

MagE 

H13, N13 

CDFS-6482 

03:32:49.34 

-27:59:52.35 

76 + 52 

15.4 + 8.1 

MMIRS (HK) 

MagE 

H13, N13 

COSMOS-08501 

10:01:16.80 

+02:05:36.26 

280 ± 30 

8.8 + 1.1 

NIRSPEC (K) 

MagE 

N13 

COSMOS-30679 

10:00:29.81 

+02:18:49.00 

87 + 7 

8.5 ±0.7 

NIRSPEC (H and K) 

MagE 

H13, N13 

COSMOS-13636 

09:59:59.38 

+02:08:38.36 

73 + 5 

11.3 + 0.5 

FMOS (H) 
NIRSPEC (K) 

MagE and LRIS 

H13, N13, S14 

COSMOS-43982 6 

09:59:54.39 

+02:26:29.96 

130 ± 12 

11.0 + 0.5 

MMIRS (HK) 

MagE and LRIS 

H13, N13, S14 

COSMOS-08357 

09:59:59.07 

+02:05:31.60 

47 + 8 

0.5 ±0.1 

FMOS (H) 

LRIS 

S14, N15 

COSMOS-12805 

10:00:15.29 

+02:08:07.50 

34 + 6 

2.6 ±0.3 

FMOS (H) 

LRIS 

S14, N15 

COSMOS-13138 

10:00:02.61 

+02:08:24.50 

40 + 10 

0.4 ±0.1 

FMOS (H) 

LRIS 

S14, N15 

COSMOS-38380 

09:59:40.94 

+02:23:04.20 

137 + 15 

2.6 ±0.3 

FMOS (H) 

LRIS 

S14, N15 

SXDS-10600 

02:17:46.09 

-06:57:05.00 

58 + 3 

1.9 ±0.1 

FMOS (H) 

LRIS 

S14, N15 

SXDS-10942 

02:17:59.54 

-06:57:25.60 

135 ± 10 

0.3 ±0.0 

FMOS (H) 

LRIS 

S14, N15 


Note. — (1) Object ID; (2), (3) Right Ascension and Declination; (4), (5) rest-frame Lyo; EW and luminosity derived from narrow- and 
broadband photometry; (6) Instruments and filters used for the NIR observations; (7) Instruments used for the optical observations; and (8) Source 
of the information 

a H13: lHashimoto et al.l (|2013D ; N13: IN aka lima et all (|2013l) : S14: iShibuva et al.l d2014bT) ; N15: Nakajima et al (2015, in preparation) 
b AGN-like object 
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is adopted as the line center, and the wavelength range 
encompassing half the maximum flux is adopted as the 
FWHM. The standard deviation of the distribution of 
measurements from the 10 3 artificial spectra is adopted 
as the error on the line center and FWHM. When mul¬ 
tiple lines are detected, we adopt a weighted mean value 
of them. A summary of the measurements are listed in 
the columns 2 and 3 of Tabic [2] All redshift (FWHM) 
values are corrected for the LSR motion (instrumental 
resolution). When the line is unresolved, the instru¬ 
mental resolution is given as an upper limit. The mean 
FWHM value for a sample of eight objects with a mea¬ 
surable velocity dispersion is FWHM(neb) = 129±55 km 
s _1 , which is smaller than that of LBGs, FWHM(neb) 
= 200 —250 km s" 1 dPettini et al.ll200HlErb et al.ll2006at 
iKulas et all 120121). This_ is consistent with the recent 
results bv lErb et all (|2014H . who have found that the 
median FWHM(neb) of 36 z ~ 2 LAEs is 127 km s _1 . 
These results indicate that LAEs have smaller dynamical 
masses than LBGs. 

3.2. Two Component [Om] Profiles 

Among the nebular emission lines we have obtained, 
while most objects show normal symmetric Gaussian pro¬ 
files, COSMOS-13138 and SXDS-10600 show an asym¬ 
metric [Oiii] profile with a secondary blueshifted and 
redshifted component, respectively (see Figure [2]). Such 
a profile has been reported in various objects: both lo¬ 
cal and high-z star- f orming galax i es an d ULIRGs (e.g., 
Shapiro et al.l 120091: IGenzel et al.l 120111 : iNewman et al.1 
2012tlSoto et al.ll2012ll a hig h-z Oxygen-Two Blob ([On] 
blobl dHarikane et al.l 1 2 01 111 . and a few Lyman-Alpha 
Blobs (LABsf (lYang et al.l 12014( 1. However in LAEs, 
there has been no study which reports its presence. 

Aforementioned studies apply a two Gaussian compo¬ 
nents fit with a narrow and broad components to the line. 
To examine the presence of two components, we also per¬ 
form a fit with two Gaussians. We have six parameters: 
fluxes, line centers, and FWHMs for both components. 
We require that the widths of both components are larger 
than the spectral resolution, and that the broad compo¬ 
nent has a larger FWHM than the narrow component. 
Best fit parameters are determined through minimum 
X 2 realizations, and the parameter range satisfying y 2 < 
Xmin + 1 is adopted as the error, where Xmin denotes the 
minimum x' 2 value. The results are listed in Table [3] For 
each object, both components are significantly detected 
with > 4 <t, demonstrating that some fraction of LAEs 
have two-component line profiles. 

The velocity offsets of the two components are 104 ± 11 
km s” 1 (COSMOS-13138) and 115 ± 8 km s" 1 (SXDS- 
10600), respectively. 

The FWHM values of the broad component after cor¬ 
rection for instrumental resolution are 70 ±50 and 80 ±30 
km s _1 . These are much smaller than those of the star 
forming galax ies at z ~ 2 (FWHM = 300 — 1000 km s _1 , 
IGenzel et al.ll20lTh . and slightly smaller than those of the 
[On] blob (FWHM = 120—130 km s -1 ) of lHarikane et alJ 
d20ll and the LA Bs (FWHM = 100 - 280 km s" 1 ) of 
lYang et all ( 2014 4. Our small values exclude the possi¬ 
bility of the broad component or igi natin g from an AGN 
activity (cf., lOstcrbrock fe Ferlandl |2006l) or a powerf ul 
outflow driven by a starburst (e.g.. iShapiro et aHl 120091 : 


IGenzel et al.ll201 ll INewman et alJl20T2l) because in these 
cases, the FWHM of the broad component should be as 
large as ~ 300 — 1000 km s -1 . It is possible that the two 
component li nes orig i nate from tw o large star-forming 
regions ( e.g.. lHa rikane e t al.ll2014fl or mergers. As dis¬ 
cussed in lHarikane et all ( 20141) . the velocity offset of the 
two components, ~ 100 km s _1 , may be due to a rotation 
of the objects. 


3.3. Lya Profile with a Blue Bump 

Wh ile the majority of Lya profiles are sing le-peaked 
(e.g., iShaplev et al .1 [20031 : iSteidel et al.l 1201011 . a frac¬ 
tion of Lya profiles are known to be multiple-peaked 
(e.g.JRauch et alll2008l: lYamada et alJl2012t IKulas et alJ 

120121 1. In particular, we shall refer to a secondary 
small peak blueward of the systemic redshift as “the 
bluebump” (see the case 2 profile in Figure 12 in 
IVerhamme et al.ll2006l) . Theoretical studies have shown 
that the blue bump is a natural outcome of the ra¬ 
diative transfer i n a lo w spee d galactic outflow (e.g., 
IZhcng fc Miralda-Escudcl 1201) 2 :■. 

We consider a blue bump to be detected if there ex¬ 
ists an excess emission blueward of the systemic redshift 
above 3er noise of the local continuum. We detect a 
blue bump of five objects; the MagE ones of CDFS-3865 
and COSMOS-43982, and the LRIS ones of COSMOS- 
12805, COSMOS-13138, COSMOS-43982, and SXDS- 
10942 (the column 4 of Table [2]). The position of the 
blue bump is designated by a blue arrow in Figure 1. 

The frequency of blue-bump objects in the sample is 
~ 40% (5/12). There are four LAEs in the literature 
that h ave a blue bump : one a mong the two LAEs stud¬ 
ied in iMcLinden et all (120111) and all three LAEs stud¬ 
ied in iChonis et al.l (120131) . For the total sample of 17 
LAEs, the frequency is calculated to be ~ 50% (9/17). 
Note that this is a lower limit due t o the l i mited spec - 
tral resolution. On the other hand, IKulas et al.l (120121 1 
have studied 18 z ~ 2 — 3 LBGs with z sys measurements 
which are preselected to have multiple-peaked Lya pro¬ 
files. They have argued that ~ 30% of the parent sample 
are multiple-peaked and that 11 out of the 18 objects 
have a blue bump, indicating that the blue bump fre¬ 
quency in LBGs is ~ 20% (~ 30% x 11/18). These re¬ 
sults imply that the blue bump feature is slightly more 
common in LAEs than in LBGs although a larger sam¬ 
ple observed at higher spectral resolution is needed for a 
definite conclusion. 


3.4. Lya Velocity Properties 

We derive three velocity offsets related to the Lya line: 
the velocity offset of the main red peak of the Lya line 
with respect to the systemic redshift, 


Av 


Lya,r 


~Lya:,r -^sys 


l 


(i) 


^sys 


that of the blue bump of the Lya line with respect to the 
systemic redshift, if any, 


Av- 


Lya,b — C - 


-^Lyoijb Z& 


1 ~h Z, 


( 2 ) 


sys 


and that of the two peaks, 

^^peak — ^^Lyc»!,r ^^Lyo!,b) 


( 3 ) 
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TABLE 2 

Summary of the observed spectroscopic properties of the sample 


Object 

•Zsys 

FWHM (neb) 

Blue Bump 

A^Lya^r 

^^Lya,b 

A^peak 

EW (Lya)spec 

Sw 



(km s -1 ) 


(km S -1 ) 

(km s _1 ) 

(km s _1 ) 

(A) 


(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

CDFS-3865 

2.17242 ±0.00016 

242 ± 31 

yes 

245 ± 36 

-352 ± 59 

597 ± 67 

40 ±2 

8.8 ±0.3 

CDFS-6482 

2.20490 ± 0.00042 

99+ 66 

yy -99 

no 

118 ±48 

_ 

_ 

26 ±2 

6.6 ± 1.7 

COSMOS-08501 

2.16161 ± 0.00042 

< 200 

no 

82 ±40 

_ 

_ 

10 ± 1 

2.2 ±2.7 

COSMOS-30679 

2.19725 ± 0.00020 

92 ±45 

no 

290 ± 33 

_ 

_ 

10 ± 1 

3.1 ± 1.0 

COSMOS-13636 (MagE) 

2.16075 ± 0.00019 

73 ± 5 

no 

146 ± 25 



23 ± 5 

5.3 ± 1.0 

COSMOS-13636 (LR.IS) 

2.16075 ± 0.00019 

73 ±5 

no 

161 ± 18 

- 

- 

26 ± 1 

6.2 ±0.5 

COSMOS-43982 (MagE) 

2.19267 ±0.00036 

325 ± 36 

yes 

117 ±53 

-297 ±57 

414 ± 78 

24 ± 17 

7.9 ± 1.3 

COSMOS-43982 (LRIS) 

2.19267 ±0.00036 

325 ± 36 

yes 

155 ± 40 

-165 ±90 

320 ± 98 

42 ±3 

-4.2 ±0.6 

COSMOS-08357 

2.18053 ± 0.00031 

< 136 

no 

106 ± 71 

_ 

_ 

19 ±3 

-4.2 ±7.9 

COSMOS-12805 

2.15887 ±0.00024 

110 ± 16 

yes 

171 ± 25 

-605 ± 114 

776 ±117 

24 ± 1 

8.9 ±0.7 

COSMOS-13138 

2.17914 ±0.00012 

63 ±6 

yes 

191 ± 59 

-214 ±87 

405 ± 105 

46 ± 11 

-1.6 ±4.8 

COSMOS-38380 

2.21245 ± 0.00015 

99 ±9 

no 

338 ± 21 

_ 

_ 

73 ± 7 

2.5 ±0.8 

SXDS-10600 

2.20922 ± 0.00014 

55 ± 28 

no 

186 ± 13 

_ 

_ 

44 ± 1 

11.7 ±0.2 

SXDS-10942 

2.19574 ±0.00025 

< 136 

yes 

135 ± 10 

-374 ±41 

556 ± 66 

94 ± 10 

1.3 ±0.3 


Note. — The symbol indicates we have no measurement. (1) Object ID; (2) Systemic redshift derived from the weighted mean of the nebular 
emission redshifts; (3) Weighted mean FWHM of nebular emission line; (4) Presence of a blue bump emission in the Lya profile; (5) Velocity offset 
of the Lya main red peak with respect to z sys ; (6) Velocity offset of the Lya blue-bump with respect to -z sys ; (7) Separation between A^Lyc^r and 
AfLycc.b; (8) Rest-frame Lya EW derived from spectroscopy; and (9) Weighted skewness of the Lya line. 




Fig. 2.— [Om] A5007 spectrum of COSMOS-13138 (left panel) and SXDS-10600 (right) taken by FMOS. The black solid lines are the 
observed spectra. There is a broad component in addition to a narrow component. The red solid lines denote the best-fit functions for the 
whole lines, and the red dashed lines represent those of the individual components. 


where z sys , ZLya,r, and ^Lya.b represent the systemic red- 
shift, the Lya redshift of the main red peak, and that of 
the blue bump, respectively. 

3.4.1. Lya Main Red Peak Velocity Offsets, AvL, ya ,i 

We estimate the AuL yQ , r value using a Monte Carlo 
technique in a similar manner to that in < 13.11 First, for 
each object, we measure the ler error in the Lya spec¬ 
trum set by the continuum level at the wavelength longer 
than 1216A. Then we create 10 3 fake spectra converted 
to velocity space by simultaneously perturbing the flux 


at each wavelength and the systemic redshift listed in Ta- 
ble[2]by their ler errors. Finally, we measure the velocity 
at the highest flux peak. The mean and the standard de¬ 
viation value of the distribution of 10 3 measurements are 
adopted as the AuL yQ , r and its error, respectively. The 
derived Ai>L y a,r values are listed in the column 5 of Table 
[21 ranging from 82 km s _1 to 338 km s _1 with a mean 
value of 174 ± 19 km s _1 . In most c ases, these values 
are co nsiste nt with those measure d in lHashimoto et ahl 
(|2013H and iShibuva et al.l (l2014bfl within ler, however, 
they are not for COSMOS-08357 and COSMOS-12805. 
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TABLE 3 

Summary of the Two-Component Fits for [Oiii] line 


Object 

/narrow 


■^narrow 


FWHM narrow 

FWHM broad 


(10 _ 17 erg s -1 cm -2 ) 

(10 -17 erg s _1 cm -2 ) 



(km s -1 ) 

(km s x ) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

COSMOS-13138 

4.7 ±0.2 

0.9 ±0.2 

2.17931 ±0.00002 

2.17821 ±0.00012 

30 ±20 

70 ±50 

SXDS-10600 

14.9 ±0.6 

4.2 ±0.6 

2.20915 ± 0.00002 

2.21038 ± 0.00008 

20 ±30 

80 ±30 


Note. — (1) Object ID; (2), (3) Fluxes of the narrow (/narrow) and the broad (/broad) components. Note that the values are not corrected for the 
slit loss; (4), (5) Redshifts of the narrow (^narrow) and the broad (.Zbroad) components; (6), (7) FWHM measurements of the narrow (FWHM nar row) 
and the broad (FWHMbroad) components. 


This is due to the fact that these studies have applied 
a symmetric/asymmetric profile fit to the Lya line. In 
Figure 1, we show the two AuLya.r values derived from 
the Monte Carlo and the profile fit technique as the or¬ 
ange and green line segments, respectively. For the sake 
of consi stency in the definition of the AuT, vr ,.r in the shell 
model (lYerhamme et all 120061 ISchaerer et all 120111) . we 
adopt here the new measurements. We note that our 
discussion is unchanged even if we adopt the previous 
AuLya.r values. 


The Au T, vry r value has been measured in more than 
60 LAEs (LMcddn den_etalJ 12Q1II: [Fin celstejnet_aL| 12Q1II; 


Hashimoto et all 20131 Guaita et a ll 20131 IChonis et *al 


2013 r iShibuva et al.l I2014bt ISong et all 12014 lErb et al 

201411 . These studies have shown that LAEs at z ~ 2 — 3 
have a mean AuLya.r of ~ 200 km s^ 1 , which is sig¬ 
nificantly smaller than that of LB Gs at a sim ila r red- 
shift, Avi, V n r ~ 400 km s -1 (e.g., iSteidel et al.1120101 
iRakic et~ahl 120111 : iKulas et~ahl 1201211 . The left panel of 


Figure [3] represents the histogram of AuL yc *,r for the 


12 LAEs (14 s pectral stud i ed in this study and 18 
LBGs given bv IKulas et al.l (12012H . We carry out the 


Kolmogorov-Smirnov (K-S) test for the two populations. 
The resultant probability is 10~ 6 , indicating that AuLy Q ,r 
is definitively different between LAEs and LBGs. 


3.4.2. Lya Blue Bump Velocity Offsets, AuL ya ,b 

For each detected blue bump in 43.31 we measure 
AuLya.b value in the same manner as for AuL y a, r - We 
obtain AuLya.b = -352 ± 59 km s” 1 (CDFS-3865), 
-297 ±57 km s” 1 (MagE-COSMOS-43982), -605 ± 114 
km s -1 (COSMOS-12805), -214±87 km s" 1 (COSMOS- 
13138), -165 ± 90 km s" 1 (LRIS-COSMOS-43982), and 
—374±41 km s -1 (SXDS-10942) as listed in the column 6 
of Table [2] We have obtained two different measurements 
for COSMOS-43982 due to the spectral resolution ef¬ 
fect, however, they are consistent with each other within 
la (see also 44.3.31) . We combine our AvL ya ,b measure¬ 
ments with those in the four aforementioned LAEs with 
a blue bump to construct a large sample of LAEs with 
a bl ue bump consisti ng of 9 objects (10 spectra): one 
from IMcLinden et al.l ()201 ll) with AuL yQ ,b = —454 km 
s -1 and three from IChonis et all (120131) with AuLya.b = 
— 127, —250, and —236 km s -1 . The mean Ai>L ya ,b value 
of the large sample is AuL ya ,b = —316 ± 45 km s -1 , 
which is consistent with that of 11 L BGs wit h a blue 
bump, AuLya.b = —367±46 km s -1 (Kulas et al.ll2012l ). 
We calculate the K-S probability to be 0.3901, indicat¬ 
ing that LAEs’ AuLya.b values are comparable to LBGs’. 
The middle panel of Figure [3] shows the A^Lya.b distri¬ 


bution for the LAE and LBG samples. 

We check if our conclusion remains unchanged even 
if the spectral resolution ef f ect is taken into account. 
The sample by IKulas et al.1 (|2012l) has been obtained 
with three settings: 300-line grating, 400-, and 600- 
line grisms, corresponding to a spectral resolution of 
R ~ 600, 800, and 1300, respectively. We compare the 
mean At>L y a,b value of our four LAEs taken by LRIS 
(.R ~ 1100) and that of six LBGs with a blue bump ob¬ 
tained at a similar resolution (R ~ 1300). The resultant 
mean AtiL ya ,b values for LAEs and LBGs are —340 ± 99 
and —356 ± 70 km s -1 , respectively, and the K-S proba¬ 
bility is 0.9238. Thus, we obtain the same conclusion. 

3.4.3. Velocity Offsets Between the Main Red Peak and the 
Blue Bump, Au pea k 

Finally, for each of the spectra with a blue bump, we 
measure the velocity offset between the red and blue 
peaks: Au pea k = 597±67 km s" 1 (CDFS-3865), 414±78 
km s' 1 (MagE-COSMOS-43982), 776 ± 117 km s" 1 
(COSMOS-12805), 405± 105 km s" 1 (COSMOS-13138), 
320 ± 98 km s” 1 (LRIS-COSMOS-43982), and 556 ± 66 
km s -1 (SXDS-10942), as listed in the column 7 of Table 
[2] In order to make a large sample with Ai> pea k mea¬ 
sured, we utilize again the four LAEs wit h the bl ue bump 
from t he literature: one LAE studied in IMcLinden et all 
(120111) with A-Un^ak = 796 km s -1 and three LAEs stud¬ 
ied in IChonis et all (I2013D with Au pea k = 300,425, and 
415 km s -1 . The mean value of the nine objects (ten 
spectra) is Au pea k = 500 ± 56 km s -1 , which is signifi¬ 
cantly smaller than the value derived for 11 LBGs with 
a blue bump, Ay nm k = 801 ± 41 km s -1 (Group I in 
IKulas et al.ll20l2 ). The K-S probability is calculated to 
be 0.00636, indicating that LAEs and LBGs have dis¬ 
tinctive Aup ea k values. See the right panel of Figure [3] 
for their distributions. 

We examine the spectral resolution effect exactly the 
same manner as in 43.4.21 The mean Au pea k value of the 
four LAEs taken by LRIS (R ~ 1100) and that of the 
six LBGs with a blue bump obtained at a similar spec¬ 
tral resolution (i? ~ 1300) are Ai> pea k = 514 ± 100 and 
778 ± 59 km s^ 1 , respectively. In conjunction with the 
K-S probability, 0.09524, we conclude that LAEs have a 
significantly smaller Au pea k value than that of LBG even 
at the same spectral resolution. Our finding is recentl y 
supported bv lHenrv et al.l (120151 ) and lYang et al.l(l2015l) . 
who have examined Lya velocity properties and their re¬ 
lations to the Lya escap e fraction for local galaxi es called 
“Green Peas” galaxies dCardamone et al.Tl2009D . They 
have found that the Lya escape fraction is higher for 
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objects with smaller Ari pea k- 

In summary, we have derived three Lya velocity off¬ 
sets, AvLy a ,r, AriLya.b; and Au pea k- While we need a 
larger sample of objects with a blue bump for a definite 
conclusion, we find that LAEs have a smaller (compa¬ 
rable) AvLyc^r (A^Lya.b) value relative to LBGs, which 
makes their Au pea k value also smaller than that of LBGs. 

3 . 5 . Other Physical Quantities 

In this section, we describe other physical quantities 
related to with this work. We describe metal absorption 
line properties in t!3.5.1l SED fitting properties in mm 
and morphological properties in i !3.5.3l 

3 . 5 . 1 . Metal Absorption Line Properties 

Low ionization state (LIS) metal absorption lines en¬ 
code information on cold neutral gas in galaxies. The 
mean blueshift of LIS absorption lines with respect to the 
systemic velocity, Au a h s , gives the average speed of the 
galactic outflow (e .g., iPettini et al.ll200il iShaplev et al.l 
120031 (Martini[20051) . In the following sections, we com¬ 
pare the Au a bs values of our LAE sample with the results 
fr om Lya radiat i on tran sfer fitting. 

iShibuva et al.l (l2014bl j have detected several LIS ab¬ 
sorption lines in a few narrowband-selected LAEs on 
the individual basis. The derived mean blue shifts are 
A Vahs = —130±70 km s -1 (COSMOS-13636), -170±50 
kms" 1 (COSMOS-128 05), and -260±60 km s" 1 (SXDS- 
10600). Additionally, lHashimoto et al.l (1201311 have de¬ 
tected several LIS absorption lines in a stacked spec¬ 
trum of four LAEs: CDFS-3865, CDFS-6482, COSMOS- 
13636, and COSMOS-30679. The mean blueshift of the 
LIS metal absorption lines is Au a b s = —102 ± 65 km s -1 . 
These values are listed in the column 2 of Tabic 2) 

3 . 5 . 2 . SED Fitting Properties 

In this study, we utilize SED fitting results of the sam¬ 
ple, in particular, stellar dust extinction, E(B — V”)*, and 
stellar mass, M*. In the following sections, we compare 
the E(B — V )* values with the results from Lya radiation 
transfer fitting, and investigate the correlation between 
the Lya profile trends and M* 

SED fitting results fo r the MagE (LRIS) o bjects 
have been presented i n lHashimoto et all (1201311 and 
iNakaiima et al.l (|2013l) (IShibuva et al.ll2014blb For the 
detai l pro cedure of the fitting, we refer the reader to 
lOno et ahl (l2010al Tbh . The derived E(B — V )* and M* 
values are listed in the columns 3 and 4 in Table [4] The 
former range from E(B — V)* = 0.04 to 0.40 with a 
mean value of E(B — V )* = 0.16, and the latter from log 
M*/Mq = 7.7 to 10.8 with a mean of log M*/Mq = 9.3, 
respectively. 

3 . 5 . 3 . Morphological Properties 

In the following sections, we use three morp h ological 
propert ies studied for z ~ 2.2 LAEs in IShibuva et al.l 
( 201 la :: the presence of a merger, the spatial offset be¬ 
tween Lya and stel lar-continuum emissio n peaks, ^Lya, 
and the ellipticity. IShibuva et alJ ( 2014al) have utilized 
/gi 4 and Hiqq data taken with ACS and WFC3 on HST 
to examine rest-frame UV and optical morphologies, re¬ 
spectively. Among the objects presented in this study, 


the rest-frame UV images of the eight COSMOS objects 
have been investigated in IShibuva et al.l (l2014aT) . 

The presence of a merger has been examined with 
two methods: the close- pair methodle.g.. iLe Fevre et aD 
l2000t lLaw et all 120121 ) and the morphological index 
method, especially C AS system (lAbraham et alJ 119961 : 
iConselice et al.ll2000f) . In IShibuva et all (l2014af) . the for- 
mer method has been applied to objects with /§ 14 < 26.5, 
which is the case for all the COSMOS objects presented 
in this study except for COSMOS-13138. The result is 
that two objects, COSMOS-13636 and COSMOS-12805, 
are mergers, while the remaining seven are not. On the 
other hand, the latter method has been done for objects 
with Jgi 4 < 25.0 and a half lig ht radi u s, r s , larger than 
0."09. The reason why IShibuva ct all (l2014al l have lim¬ 
ited the sample for the latter method is to obtain re¬ 
liable values of the indices. This is the case for three 
COSMOS objects presented in this study, COSMOS- 
13636, COSMOS-43982, and COSMOS-38380. The re¬ 
sult is that none of the three is a merger. The two re¬ 
sults for COSMOS-13636 are not consistent with each 
other because we have used two different methods. Thus, 
among the eight COSMOS objects, COSMOS-13636 and 
COSMOS-12805 may be a merger (the column 5 of Table 

13). 

The Lya spatial offset, (5Lyaj has been examined by 
performing source detections with SExtractor for Sub¬ 
aru NB387 and HST Jgi 4 images. While compact ob¬ 
jects with a symmetric UV light profile tend to have a 
small <5Lya value, objects with an asymmetric, disturbed 
UV light profile likely to have a larg, e Si, ya value (e.g., 
iJiang ct ahll20131 : IShibuva et al.ll2014al) . Thus, this quan¬ 
tity could be a useful tracer of the Hi gas stability around 
the galaxy. The value is reliably obtained for the objects 
with Igi 4 < 26.5 and NB387 < 24.5, where the typical 
positional error in Igi 4 (NB387) is less than 0."02 (0."3). 
For the eight COSMOS objects in this study, none has a 
significant Lya spatial offset larger than the typical error 
of the <5Lya> ~ 0."36. 

The ellipticity, e = 1 — a/ 6 , where a and b are the major 
and minor axe s, is a usefu l in dicator of the galactic disk 
inclination. In lShibuva et alJ (l 2014a|), this has be en mea¬ 
sured using GALFIT software (iPeng et al.ll2003l for the 
objects with J 8 i 4 < 25.0 and r e larger than the typical 
PSF size. The former criterion, corresponding to S/N = 
30 detectio n, is needed for the reliable ellipt i city m easure- 
ments (e.g. lMosleh et al.l 120121 : lOno et al.l 1201 3l i. Only 
three objects, COSMOS-30679, COSMOS-38380, and 
COSMOS-43982, satisfy these criteria. The resultant 
ellipticity values are e = 0.24 (COSMOS-30679), 0.34 
(COSMOS-38380), and 0.49 (AGN-COSMOS-43982), re¬ 
spectively (the column 6 of Table |4j. 

4. LYo RADIATIVE TRANSFER MODEL AND FITTING 
PROCEDURE 

4.1. A Library of Synthetic Spectra 

The library of synthetic Ly a spectra used in thi s 
study has been described in iSchaerer et alJ ( 2011 ). 
Lya radiation transfer has been computed with McLya 
(IVerhannne et all 120061) through spherically symmetric 
expanding shells of homogeneous and isothermal neutral 
hydrogen gas. The shell is describe by four parameters: 

• the radial expansion velocity, T4xp> 
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Av L y a r [kms '] Av Lycx ,b [kms -1 ] Av peak [kms 1 ] 

Fig. 3. — Histograms of A^Lyoi.r (left panel), At»L yQ , b (middle), and Ai;p ea k (right) for the LAEs studied in this paper and the literature 
(|McLinden et alll2QlU : FChonis et al.|[2013l ) and LBGs given bv IKulas et al.l (|20121 ). 


TABLE 4 

Summary of the Other Physical Properties 


Object 

^^abs 

log (M,) 

E{B - V), 

merger 

€ 





(pair, CAS) 



(km s -1 ) 

(M 0 ) 

(mag.) 



(i) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

CDFS-3865 

(-102 ±65) 

9-50to.M 

0 14+ 0 - 00 

U,J - - 0.00 

— 5 “ 

- 

CDFS-6482 

(-102 ±65) 

9.80^;- 

0-15^-.O2 

“j “ 

- 

COSMOS-08501 

- 

7.84+J;^ 

0.08+qqI 

no, - 

- 

COSMOS-30679 

(-102 ±65) 

' -0.52 

0 24 + ^ 04 

u>z ^-0.04 

no, - 

0.24 

COSMOS-13636 

-130 ±70 (-102 ±65) 

9.12+^ 

n i c+b-ui 

U-Io 0.01 

yes, no 

- 

COSMOS-43982 

- 

10.80^QQg 

0 40 , ^ uu 

u.4tu_ 0> oo 

no, no 

0.49 

COSMOS-08357 

- 

Q 91 

y - Zi -0.40 

0 i 4 + u - u£) 

u.14 o 05 

no, - 

- 

COSMOS-12805 

-170 ±50 

y.^o 17 

0 i «+ u - u2 

u - iO - 0.02 

yes, - 

- 

COSMOS-13138 

- 

Q 4Q+0.22 

y-4o 0.20 

0 i 9+^- U4 

U ‘ iy -0.04 

- 

- 

COSMOSO-38380 

- 

10.06+^5 

0 1 Q+U.U 2 

u.lo 0 01 

no, no 

0.34 

SXDS-10600 

-260 ± 60 

9 4U -|u - ut) 

0.04 

0 f)5+ u - uu 

U.UD_o 01 

- j - 

- 

SXDS-10942 

- 

7 yo+U-LL 

' • /o -0.08 

U -U4±_0 02 

— 5 “ 

- 


Note. — The symbol indicates no measurement. (1) Object ID; (2) Mean velocity offset of LIS absorption lines with respect to -z sys ; (3) 
Stellar mass estimated from SED fit ting; (4) Stella r dust e xtinction estimated from SED fitting; (5) Presence of merger examined via close-pair 
method and CAS system studied in lShibuya et all i2014aD ; (6) Ellipticity defined as e = 1 — b/a , where a and b are the major and minor axes, 
respectively. 


• the neutral hydrogen column density along the line 
of sight, N m , 

• the Doppler parameter, b 1 describing the thermal 
and turbulent motion in the shell, 

• and the dust absorption optical depth at the Lya 
wavelength, r a , related to the gas dust extinction 
by E(B — V) gas ~ (0.06...0.11)r a , where the lower 
and higher values in the parenthesis correspond to 
the a ttenuation law for starbursts (ICalzetti et al. 
200C) and the Galactic extinction law dSeaton 
197911 . respectively. 


For a comparison with the observed data, each rest-frame 
model has been shifted using the systemic redshift 2 sys 
values listed in Table [2] To reflect the 2 sys uncertainty, 
we have allowed the observed Lya spectra to shift rel¬ 
ative to the velocity zero point within the error. Thus, 
combinations of seven free parameters are fitted to the 
data. 

This library of Lya spectra has been successfully 
used to reproduce various observed Lya line profiles 
of z > 3 LBGs, from strong emission to broad ab¬ 
sorpt i on llVerham me et all 120081 iSchaerer fe Verhamme 
2008 : iDessauees-Zavadskv et al.M201d : IVanzella et al. 
201C ; iLidman et al.ll2012lf . 


The Lya source is located at the center of the shell. The 
intrinsic (i.e., before being affected by the radiative trans¬ 
fer effect) spectrum is a Gaussian Lya line plus a flat 
continuum, and is characterized by two parameters : 

• the Lya equivalent width, EWi n t(Lya), 

• and the full width at half maximum, 
FWHM int (Lya). 


4.2. Fitting of Observed Spectra 

To perform a statistical comparison between the ob¬ 
served and modeled Lya line profiles, we calculate the 
X 2 values for each of the possi ble combinations of the pa¬ 
rameters for each galaxy (cb. lGhonis et al.ll201.3 i. Since 
model spectra are normalized and at an infinite spectral 
resolution, two steps are needed before the % 2 calcula¬ 
tion. First, we normalize the observed spectra using the 
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continuum level estimated at wavelengths longer than 
1216A. Second, each model Lya spectrum has been con¬ 
volved with a Gaussian whose FWHM corresponds to 
spectral resolutions: 

FWHM = c/R, (4) 

where c is the speed of light. 

We note that our fitting technique gives exactly the 
same statistical weight to all data points of the contin¬ 
uum and the Lya line. Finally for the sake of consistency, 
for each object we calculate y 2 in the wavelength range 
from —3 x FWHM Q bs(Lya) to +3 x FWHM Q b s (Lya) 
around the Lya line center. 

In Figure HI we demonstrate how the best fit, and its 
associated errors, are found using y 2 values. To do this, 
examples of the fit to 14 xp are shown for well and poorly 
constrained objects. In the left panels of this figure, one 
can see a broad range of V exp values with low reduced 
y 2 for COSMOS-08357 whose Lya S/N ratio is ~ 11, in 
comparison to CDFS-3865 with a Lya S/N of ~ 98. To 
measure median and ler values, we convert y 2 values into 
probabilities using the formula, p oc exp(—y 2 /2) for each 
five 2D parameter set (I4xp vs. Ahi, I4 xp vs. r a , T4 xp vs. 
6, V/xp vs. FWHMi n t(Lya), and V exp vs. EW int (Lya)). 
After normalizing them so that the total probability is 
unity, we draw a probability (PDF) and a cumulative 
density function (CDF) as shown in the middle and the 
right panels, respectively. Finally, we adopt the values 
where the CDF value satisfying CDF = 0.50,0.16, and 
0.84 as the median and ±1 <t, respectively. Performing 
this for each five 2D parameter set results in five median 
and ±1 <t values. As can be seen, all the five median and 
±1 a values are consistent with each other for CDFS- 
3865, whereas those are not for COSMOS-08357. In the 
latter case, we adopt the average of the five median and 
±1<T values. 

4.3. Results 

We show the reproduced Lya profiles fi l4.3.1D . describe 
the derived parameters 1 EI4.3.2II . and examine the influ¬ 
ence of spectral resolution on the results ( 114.3.31) . 

4.3.1. Fitted Profiles 

Figure [5] shows the best fit model spectra with the 
observed ones. All the Lya profiles are quite well re¬ 
produced by the m odel, which seems to d iffer from the 
previo us studies bv lKulas et al.1 (|2012f) and lChonis et al.1 
(I2013T) . These authors have had difficulty reproducing 
their Lya profiles, especially the position and the flux of 
the blue bump. This might be due to model differences. 
These two studies have utili zed the uni form expand¬ 
ing sh ell m odel constructed bv I Zheng &; Miralda-Escudl 
(l2002f) and iKollmeier et al.l (1201011 . There are three ma¬ 
jor d ifferences between the models (c.f., iChonis et al.l 
120131) . First, in addition to the three common param¬ 
eters, I4x P , IVhi, and b, the model used in this study also 
includes an additional one for dust absorption. Second, 
the grid points and the phy sical range o f parameters are 
different. The model bvIZheng fe Miralda-Escudel (|2002f) 
and IKollmeier et al.l (l201(j| f has four values for each pa¬ 
rameter: V exp = 50, 100, 200, 300 km s -1 , log (Ahi) 
= 17, 18, 19, 20.3 cm~ 2 , and b = 20, 40, 80, 120 km 
s -1 , whereas the model used in this study has 12 Kx P , 


13 Ahi, and 5 b values spanning wider physical ranges. 
Finally, the intrinsic spectrum of the previous models is 
a monochromatic Lya line, while we model a Gaussian 
Lya plus a continuum. As we show in H4.3.2I and later 
sections, we infer that the key to better reproducing the 
blue bump is to assume the Lya profile to be a (broad) 
Gaussian. 

4 . 3 . 2 . Derived Parameters 

The best fit parameters are summarized in Table [5j 
We describe the mean values of the derived parame¬ 
ters, and systematically compar e them with those of 
L BGs modeled by the same c o de (iVerhannne et al.H2008l: 
Schaerer fe Verhammel 120081 : iDessauees- Zavadsky et al.1 
2010T) . For the parameter FWHMi nt (Lya), we examine 
the mean values of two subsamples, objects with a blue 
bump and those without. We have checked that there is 
no significant difference between the two subsamples in 
the other parameters. 

The mean I4 xp value of the LAEs is 148 ± 14 km s _1 , 
which is comparable to that of LBGs, 131 ± 25 km s _1 . 
This strongly disfavors the hypothesis that the small 
AuLy a of LAEs is due to their large outflow velocity. 

The most interesting parameter, Ahi, ranges from 
log(Am) = 16.0 to 19.7 cm- 2 , with a mean value of 
18.4 ± 0.4 cm -2 , which is more than one order of mag¬ 
nitude smaller than the typical log(Am) value of LBGs, 
19.8 ±0.2 (cm" 2 ). 

The mean values of r a and b are 0.9 ± 0.2 and 37 ± 10 
km s, respectively, both of which are comparable to 
those of LBGs, 0.8 ih 0.1 and 28 ± 5 km s _1 . 

FWHMi n t(Lya) values range from FWHMi nt (Lya) = 
50 to 847 km s _1 . The mean values for the whole sam¬ 
ple, the non blue bump sample, and the blue bump 
sample, are 354, 169, and 602 km s _1 , respectively. 
This shows that the blue bump objects have significantly 
larger FWHMi nt (Lya) than that found i n the non b lue 
bump objects. This trend is similar to IVerhamme et all 
(12008T) : They have found that most LBGs with a sin¬ 
gle peaked Lya profile are best fitted with moderate 
values of FWHMi nt (Lya), ~ 200 km s _1 , whereas the 
best fit FWHMi n t (Lya) values for two LBGs with a blue 
bump are greater than 500 km s _1 . These results sup¬ 
port our claim that large FWHMi nt (Lya) helps fitting 
the blue bump. We investigate if there are any obser¬ 
vational trends for the blue bump objects, and discuss 
possible mechanisms for the blue bump objects to have 
large FWHM; nt (Lya) in i l5.ll 

Since starburst activities that produce Lya photons 
should be similar between LAEs and LBGs, we expect 
comparable EW; n t (Lya) values for these two galaxy pop¬ 
ulations. The result is that the mean EWi n t(Lya) value 
of LAEs, 65 ± 18 A, is somewhat smaller than that of 
LBGs, 107 ±25 A. 

In summary, the model parameter Ahi derived in LAEs 
is more than one order of magnitude smaller than that of 
LBGs, whereas the remaining parameters are consistent 
within ler between LAEs and LBGs. 

4 . 3 . 3 . Influence of Spectral Resolution on the Fitting 
Procedure 

To investigate the influence of spectral resolution on 
the fitting results, we compare the best fit parameters 
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of the two objects observed with the two spectrographs, 
COSMOS-13636 and COSMOS-43982. As can be seen 
in Table [U the two fitting results of COSMOS-43982 are 
consistent with each other, whereas those of COSMOS- 
13636 are not, possibly owing to the large difference in 
the best-fit reduced y 2 , 1.1 and 6.2. 

Taking a closer look into these two fits, we see that the 
extremely small 1 a noise in the flux of LRIS-COSMOS- 
13636 could be a key reason for its high x 2 value. On 
the other hand, the modeled spectrum seems to be over¬ 
smoothed, leading us to infer its Lya line resolution is 
under-estimated. Indeed, it is known that the spectral 
resolution for a given line can be higher than the canoni¬ 
cal value. A combination of these factors would naturally 
cause the large resultant y 2 value, and the discrepancy 
between the different best-fit parameters at two resolu¬ 
tions. 

4.4. Degeneracy among Parameters 

In this subsection, we investigate degeneracies among 
the model parameters to understand how they affect our 
determination of the best fit parameters. First we de¬ 
scribe possible degeneracies and then statistically exam¬ 
ine them using 2D y 2 values. 

It is possible that parameters r and EW; n t(Lya) are 
degenerated as an observed profile can be reproduced 
equivalently well either assuming a weak intrinsic line 
with low dust extinction, or a strong intrinsic line with 
high dust extinction. There would also be a degeneracy 
between b and FWHMi nt (Lya) in the sense that both 
broaden the line profile. Furthermore, when there is a 
blue bump in the profile, we need either a high b or a low 
14x P to reproduce it. 

Figures 13 - 15 in the Appendix are 2D parameter grid 
maps for CDFS-3865 with the grey dots showing the en¬ 
tire grids. We use these maps and y 2 values to examine 
the actual degeneracies among the parameters. If there 
is a degeneracy between two parameters, the y 2 contour 
would be tilted and elongated. The blue grids in these 
figures show those satisfying Ay 2 <6.17 above the raw 
minimum y 2 designated by the whi te dots, i .e., the 3 
a uncertainty in the parameter set dPress et al.ill992 1. 
Thanks to the number of data points given by high spec¬ 
tral resolutions, and the relatively coarse grids, even the 
3 a uncertainty is converged into one grid. This indicates 
that there is no degeneracy that affects our determina¬ 
tion of the best fit. We have checked that this is also 
true for the rest of the sample in this study. Thus, we 
conclude that the systematic uncertainties among the pa¬ 
rameters due to the degeneracies are small, and thus do 
not affect our discussions. 

4.5. Comparison between Observation and Model 

In order to examine if the best fit parameters are rea¬ 
sonable, we compare the derived parameters with the 
observables. 

4.5.1. | AUabs | VS. Vexp 

As stated in H3.5.1I several LIS absorption lines have 
been detected in individual spectra of C OSMOS-12805, 
COSM OS-13636. and SXDS-10600 (IShibuva et all 
l2014bl) . and in a stacked spectrum of four LAEs, CDFS- 
3865, CDFS-6482, COSMO S-13636, and COSMOS- 
30679 (lHashimoto et all I2013T ). The measured blueshift 


of LIS absorption lines with respect to the systemic, 
Au a bs, is listed in Table [U Figure [5] shows a comparison 
between |Au a b s | and the best-fit expansion velocity, 
14xp- For the stacked spectrum, we plot the mean !4xp 
value of the four LAEs, 163 ±25 km s -1 . While there are 
only four data points, |Au a b s | and V exp are in excellent 
agreement with each other. 

4.5.2. E(B-V), vs. r a 

The stellar dust extinction values, E(B — E)*, 
for the sample have been derived in previous 
studies (lHashim oto et al.l 120131; iNakaiima et al.1 120131 
IShibuva et al.l 12014hf l Isee H3.5.1|) . Figure [7] compares 
them with gas dust extinction, E(B — H) gas , derived as¬ 
suming the relation: 

£(H - F) gas « 0.10r a . (5) 

Dotted and dashed lines corresp ond to empirical rela¬ 
tions E(B — V)„ = E(B — (lErb et al.ll2006alf and 

E(B - E)* = 0.44 E{B - C) gas flCalzetti et af .ll20 0 0i) , re- 
spectively, for host galaxies. As iKashino et al.l ~ (l2013f ) 
have shown, the difference between E(B — E)* and 
E(B — E) gas becomes smaller for higher -z galaxies. 

In this study, we expect that data points are located 
below these relations. This is because E(B — H) gas ob¬ 
tained from Lya modeling is gas dust extinction for out¬ 
flowing shells , which should be smaller than that for 
host galaxies. The figure shows that half of the sample 
roughly lie between the two lines, while the rest of the 
sample show low E(B — V) ^ values. A simil a r tren d has 
been found in Figure 12 of lVerhamme et all (120081 ) who 
have compared E(B — V )* and E(B — H) gas for z ~ 3 
LBGs. They have assumed two different star formation 
histories (SFHs) in deriving E{B — V )*: a constant SFH 
indicated by red triangles and an exponentially decreas¬ 
ing SFH indicated by blue open circles, the former of 
which is the same as that assu med in this study . Both 
our data and the red triangles in lVerhamme et ahl (l2006l l 
are similarly distributed in the sense that half of the sam¬ 
ple has comparable extinction values and the rest has low 
E(B — H) gas values. 

4.5.3. FWHM(neb) vs. FWHM int (Lya) 

Figure [8] plots the observed FWHM of nebular emis¬ 
sion lines, FWHM(neb), versus modeled FWHM of 
the intrinsic (i.e., before being affected by the radia¬ 
tive transfer effect) Lya line, FWHM; nt (Lya). Assum¬ 
ing that both Lya and nebular emission lines originate 
from Hu regions, the two FWHMs should be similar. 
However, FWHMj nt (Lya) is systematically larger than 
FWHM(neb). Additional scattering of Lya photons in 
an Hn region due to residual Hi atoms in it may be at 
work. Assuming a static Hn region with a neutral hydro¬ 
gen column density of log(Nni) < 17.0 cm -2 , correspond¬ 
ing to an unity optical depth f or ionizing photons, r lon < 
1 (cf., IVerhamme et al.l 12015 :). FWHM int (Lya) can be 
broadened by 200 km s -1 compared to FWHM(neb). As 
can be seen from Figured while this additional broaden¬ 
ing would help explain the discrepancy for the non blue 
bump objects, it is still not enough for the blue bump 
objects. We discuss some interpretations fo r th e huge 
FWHMjnt (Lya) in the blue bump objects in H5.ll 

We also perform Lya profile fitting of the blue bump 
objects with fixing FWHMj nt (Lya) = FWHM(neb). As 
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TABLE 5 

Summary of the Lyo fitting for the sample 


Object 

X? ed 

kexp 

log(JVni) 

T a 

b 

FWHM(Lya)j nt . 

EW(Lya) int . 



(km s -1 ) 

(cm -2 ) 


(km s -1 ) 

(km s -1 ) 

(A) 

(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

CDFS-3865 

3.1 

120±l\ 

19-518:1 

n n+ 0.0 

u - u -o.o 

15l“ 

846l^ 6 

35l? 

CDFS-6482 

1.3 

177±il 

19-218:1 

0 

to 

1 + 

0 0 

b b 

00 tp*. 

iolg 

o' 7 -i +38 
z ‘ 4-29 

28l? 

COSMOS-08501 

1.3 

167i 2 i*® 

18.715;® 

1 C 0 + 1-54 
l.OD—i 07 

13l“ 

2521118 

14l? 

COSMOS-30679 

1.0 

1271JJ 

19 . 518:8 

1 4Q+1-03 
1 * 4 l * 3 -0.53 

29l® 

5018® 

39lli 

COSMOS-13636 (MagE) 

1.1 

2261“ 

i 6 .ol 8:8 

o.ialgis 

12111? 

256118 1 

28l? 

COSMOS-13636 (LIRS) 

6.2 

1271“ 

18-818:1 

o.osi 8;88 

30l® 

12711® 

28l? 

COSMOS-43982 (MagE) 

1.0 

14111? 

18.2+°® 

1 1 r: + l-54 

-i.-J-o_o.89 

1211 1 

544118“ 

281k 

COSMOS-43982 (LRIS) 

1.4 

1381®® 

18-llSI 

0 02”*” 0,22 
U.UZ-0 02 

13l® 

621188 

42l? 

COSMOS-08357 

1.3 

I70I42 

i9.7i8i 

2 24”*~ 1,25 
z - Z 4 -0.95 

19l“ 

74l|2 

85118 

COSMOS-12805 

3.1 

H7l 21 

19.2I8I 

1 - 70 + 0.71 
x * 0.38 

101“ 

645188 

4218 

COSMOS-13138 

1.5 

oi+481 

Z1 -21 

16-618:? 

1 13"*” 1 ' 45 
X-XO —0.89 

15l“ 

501 II 44 

64lH 

COSMOS-38380 

2.1 

1271“ 

19-718:1 

0 69”*” 0 ' 28 
u.oy_ 0 32 

60111 

9918 

2761“ 

SXDS-10600 

6.2 

2261“ 

i 6 -ol 8:8 

1 74+ 0 - 20 

I 21 I 87 

223118 

H3l? 

SXDS-10942 

1.6 

13111s 

i6.oi8;8 

0 19 - *" 0,08 
u,xz -0.08 

60111 

4531®? 

85l“ 


Note. — (1) Object ID; (2) Reduced x 2 value of the fitting calculated as Xred = X 2 /(N — M ), where N and M denote the number of data points 
and the degree of freedom, respectively; (3) — ( 8 ) Best fit values of the radial expansion velocity, the column density of the neutral Hydrogen, the 
dust absorption optical depth, the Doppler parameter, the intrinsic Lyo FWHM, and the intrinsic Lyo EW, respectively. 


shown in Figure[9l the blue bumps are poorly reproduced 
compared to the fitting without fixing FWHMi nt (Lya). 
We examine if the derived best-fit model parameters dif¬ 
fer between the free and fixed FWHMi nt (Lya) cases. 
While there is no systematic difference for 14x P and Nm, 
we find that b (r a ) becomes large (small) in the fixed 
FWHMi n t (Lya) case. This would be related to th e in¬ 
trinsic degeneracy between them discussed in t )4.4l 

4.5.4. EW(Lya) vs. EW in t(Lya) 

Figure [12] plots the observed EW(Lya) against the 
best fit intrinsic EW(Lya) obtained from the Lya fit¬ 
ting, EWi n t(Lya). Since we have modeled Lya emis¬ 
sion lines that fall in the slit, we use EW(Lya) val¬ 
ues measured from spectra as the observed EW(Lya). 
All the data points are expected to lie above the one- 
to-one relation, EWi„t(Lya) > EW(Lya) spe c- This is 
because we have used the uniform shell model which 
does not boost EW(Lyq) unlike clumpy shell models 

fcf..lNeu feldlll99l MLaursen et alJl20l3t iDuval et al.ll20l4 

IGronke fc Diikstral 12014th As can be seen, all the data 
points satisfy the expectation within the 1 <t uncertainty. 

5. DISCUSSION 

5.1. Mystery of the Blue Bump Objects 

As described in previous sections, FWHMi nt (Lya) > 
FWHM (neb) is required to well reproduce the Lya pro¬ 
files with the blue bump objects. As seen in Figured the 


position and flux of the blue bump are poorly reproduced 
if we fix FWHMint (Lya) = FWHM(neb). 

In this section, we first examine if there are any char¬ 
acteristic properties for the blue-bump objects, and dis¬ 
cuss the origin of the large discrepancy between the two 
FWHMs. 

5.1.1. Any Difference in Properties between the Blue Bump 
and the Non Blue Bump Objects ? 

In 1)4.3.21 we have argued that, among the model pa¬ 
rameters, only FWHMint (Lya) is significantly different 
between the blue-bump objects the and non blue-bump 
objects. Here we examine the difference in stellar mass, 
Lya luminosity, morphological ellipticity, and the merger 
fraction between the two samples. 

First, it is possible that the non blue bump objects have 
faint Lya luminosities and/or small stellar masses so that 
the blue bump can not be observed. The Lya luminosity 
of the blue bump sample ranges from L(Lya) = 0.3 to 
29.8 x 10 42 erg s _1 with a mean value of 8.8 ± 5.6 x 
10 42 erg s _1 , whereas that of the non-blue bump sample 
ranges from L(Lya) = 0.5 to 15.4 x 10 42 erg s _1 with a 
mean value of 7.0 ± 2.1 x 10 42 erg s _1 . This indicates 
that the two subsamples have similar Lya luminosities. 
Likewise, the stellar mass of the blue bump sample ranges 
from log(M*/MQ) = 7.73 to 10.80, with a mean value of 
9.4 ± 0.5, whereas that of the non-blue bump sample 
ranges from log(M*/M 0 ) = 7.84 to 10.06, with a mean 
value of 9.3 ± 0.3. Thus, this possibility is unlikely. 






















14 


Hashimoto et al. 


T3 


CM 






Fig. 4.— Examples of reduced \ 2 values (left panels), converted probability density functions (PDF) (middle panels), and cumulative 
density functions (CDF) (right panels), for the parameter I4xp- The upper (lower) panels are for CDFS-3865 (COSMOS-08357). 


Second, objects with a blue bump may be more likely 
to be seen ed ge-on than those without. Recent theo¬ 
retica l studies rtVerhamme et, alJl2012t iZheng fc Wallacel 
120141) have investigated the inclination effects to the Lyct 
emissivity and profile. These studies have shown that the 
blue bump flux relative to the total Lya flux is enhanced 
with an increasing cllipticity. Indeed, Lya profiles seen 
edge-on in these simulations resemble those produced by 
the static case of the spherical shell model. This is be¬ 
cause outflowing gas is more likely to be blown out per¬ 
pendicular to the galaxy disk, reducing the relative out¬ 
flow velocity in the plane of the disk. As seen in Table 
IH there are three objects whose ellipticity has been mea¬ 
sured. Due to the small number of objects, we cannot 
determine if there is any difference between the two sub¬ 
samples. 

Finally, a, s disc ussed in iKulas et al.l (I2012D and 
IChonis et al.1 (I2013D . galaxy merging can be the origin 
of the blue bump. In this case, the redder and bluer Lya 


emission compone nts correspond to the two objects, re - 
spectively (see also lCooke et al.ll2010l ; lRauch et al.l201lD . 
However, as described in Tf.5.31 (Table 4), the merger 
fraction in our sample is quite low. 

We note here the observationa l results of lErb et ahl 
(|2010l) and iHeckman et al.l (1201 ll) . These studies have 
found that objects with a blue bump tend to have a low 
covering fraction of the neutral gas measured by LIS ab¬ 
sorption lines. Unfortunately, we cannot test this trend 
with our sample because of a too small number of objects 
with detection of LIS absorption lines. 

We conclude that there is no significant difference in 
Lya luminosity, stellar mass, morphological ellipticity, 
or the merger fraction between the two samples. A large 
sample, whose Lya and absorption line velocity proper¬ 
ties as well as morphological and stellar population prop¬ 
erties are simultaneously available, is needed to under¬ 
stand the origin of blue bumps. 
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Fig. 5. — The upper two rows of panels show the reproduced Lya line profiles (red) on top of the observed ones (grey) for the MagE 
objects, while the lower two rows of panes are those for the LRIS objects. The gray region denotes the lcr range of the observed spectrum. 
All spectra are scaled to the wavelength range from —1500 to +1500 km s -1 . 
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lAVabsl NS 1 ] 

Fig. 6 .— Vexp plotted against |Av a b s |. The diamonds are the 
three LRIS objects w hose LIS absorption li nes are detected in the 
individual spectrum (|Shibuva et al.M2014bl) . Blue color shows an 
object with a blue bump in the Lya profile, while black color de¬ 
notes those without. The magenta circle is the stacked spectrum of 
the four MagE LAEs, CDFS-3865, CDFS-64 82, COSMOS-13636, 
and COSMOS-30679 (|Hashimoto et aLll2O130 . 



Fig. 8 . — FWHMi nt (Lya) plotted against FWHM(neb) for all 
the sample. The meaning of the symbols and the colors is the same 
as in Figure|7| The solid line shows the one-to-one relation between 
the two FWHMs, while the dotted line is a relation between the 
two after taking into account an additional scattering of Lya pho¬ 
tons by residual Hi atoms in the Hu region, FWHMi n t(Lya) = 
FWHM(neb) + 200. 



E(B-V). 


Fig. 7. — E(B — V)gas plotted against E{B — V)* for all 
the sample. The circles show the MagE objects, while the dia¬ 
monds are the LRIS objects. Blue color shows objects with a blue 
bump in the Lya profile, while black color denotes those with¬ 
out. The do tted and dashed lines correspond to E(B — I/)* = 
E(B — yjo-aa dErb ct alj2006bl ) and E(B — V)* = 0A4E(B — V) Ba , s 
ICalzetti et alJl200(11 . respectively. 

5.1.2. A Possible Explanation for Large FWHMi nt (Lya) in 
Blue Bunp Objects 

In this subsection, we explore a possible explanation 
of the large discrepancy between FWHMj nt (Lya) and 
FWHM(neb) in the blue bump objects. 

It is possible that observed Lya photons are pro¬ 


duced not only from recombination of hydrogen gas 
ionized i n Hu regions, but als o from e.g., shock 
heating (lOtf-Floranes et ahl I2012D . fluorescence (e.g., 
Cantalupo et al.l 20121 20141), and/or gravitational cool- 
ing (e.g.. iDiikstra et al.ll2006D . If these are taken into 
account, the huge FWHM; nt (Lya) in the blue bump ob¬ 
jects could be explained as follows. 

Fluorescence caused by a QSO would ionize the outer- 
layer of the ISM of galaxies, and produce a large 
FWHMj n t(Lya). However, there are no QSOs around 
any of our objects. 

Gravitational cooling is another mechanism that pro¬ 
duces Lya photons. When gas inflows into the grav¬ 
itational potential well of a galaxy, the gravitational 
binding energy is converted into the thermal energy, 
which is in turn released as Lya photons (e.g., iDiikstral 
120141) . Since it occurs in both the inner and outer 
regions of the galaxy, gravitational cooling can give a 
large FWHM; nt (Lya). Furthermore, gravitational cool¬ 
ing can reproduce not only the observ ed enhanced Lya 
blue bump flux (e.g., IDiikstra et al.l I2006D. but also 
the s patially extended Lya source ( Rosdahl fe Blaizotl 
IMD, i- e -! diffuse Ly a haloes which are common fea¬ 
tures around galaxies (IS teidel et alJl201l UMatsuda et al.l 
I2012t IHaves et al .11201.11 IMomose et al.ll2014D . We note 
here that there exists a large uncertainty in model¬ 
ing the Lya emission from gravitational cooling due 
to it s difficulty and assumed observation sensitivity 
Ccf.. iFaucher-Giguere et ah I [201(1: IGoerdt et al .1 [2(71(1: 
IRosdahl fe Blaizotll2012l : lYaiima et al.ll2012L 120151) . Our 
results as well as observational results quoted above can 
be useful for future modeling. 

We conclude that not only the large discrepancy be¬ 
tween the observed FWHM(neb) and FWHM; nt (Lya), 
but also the presence of a blue bump can be simultane- 
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Fig. 9.— Ly a fittings with fixed FWHM; n t(Lya) = FWHM(neb) for the objects with a blue bump. The reproduced Lya lines (blue) 
are overlaid on the observed ones (grey). For comparison, we also plot the reproduced profiles without fixing FWHM; n t(Lya) (red). 



Fig. 10. — EWint(Lya) plotted against EW(Lya) obtained from 
spectroscopy. The meaning of the symbols and colors is the same 
as in Figure[7] Since we have assumed a uniform shell model which 
does not cause a EW(Lyo) boost, the data points are expected to 
lie above the one-to-one relation. 

ously explained if we introduce additional Lya photons 
produced by gravitational cooling. 

5.2. Origin of Small AuL ya r in LAEs 


As described in 43.41 the mean AtiLya.r of LAEs, 
~ 200 km s _1 , is significantly smaller than that of LBGs , 


Apt, v „ r ~ 400 km s 1 (LBGs: e.g.. ISteidel et all 

2010; 

Rakic et al. 20111 Kulas et al.H2012HSchenker et all 

2013, 

LAEs: e.g.. McLinden et al. 2011: 

Hashimoto et al. 

2013 

ilChonis et alJl2013klShibuva et al.l 

2014btlErb et al. 

2014 

). We have also demonstrated that some LAEs 


have an extremely sm all Atii, vr ,. r of 5, 1 00 km s -1 . 
lHashimoto et al.1 (|2013l ) and lShibuva et all (|2014bf ) have 
also found that AwLya.r correlates with SFR, velocity 
dispe rsion, ste lla r ma ss, specific SFR, and dust extinc¬ 
tion. lErb et ahl (I2014D have also found that AuLya.r cor¬ 
relates with velocity dispersion. In addition, they find 
that objects with a small AuL ya ,r have a large fraction of 
emission blueward of the systemic velocity, while the red 
wing of the Lya profile and the outflow velocity traced 
by absorp tion lines remai n unchanged. Following these 
findings. lErb et ahl (I2014D have argued that the small 
AdLya.r hr LAEs is consistent with a scenario where the 
opacity to Lya photons is reduced by a bulk motion 
and/or coverin g fraction of the ga s near the systemic ve¬ 
locity (see also lSteidel et al.lf2010f> . These results suggest 
that AiiLya.r is closely related with the physical size of 
the galaxy system. However, there are still no definitive 
conclusions why LAEs have small AuLyo.r- 
In this subsection, we explore the origin of the small 
AuLya.r in LAEs using the largest sample of LAEs 
whose high-quality spectroscopy data and several prop¬ 
erties have been obtained. There are several hy¬ 
potheses which give AuL ya ,r as small as 0 — 200 km 
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s” 1 : a unifor m shell ISM with a h igh-speed galactic 
outflow (e.g., IVerhamme et al.l I2006D . a uniform shell 
ISM with a low neutral hyd rogen column density (e.g., 
IVerhamme et al.lf2006l . l2015t) . and other m odels suc h as a 
clumpy ISM with a low covering factor f r (IHanse n fe Ohl 
120061 IDiikstra fc Kramc^l2012t lLaursen et al.l 120131 ). or 
shell models with holes/cavities, i.e . , CF < 1 (e.g., 
IBehrens et al.l 12014 IVerhamme et al.l 1201511 . We quan¬ 
titatively discuss the first two hypotheses based on our 
detailed comparison of data with uniform shell models 
1 95.2. H and 95.2.2|l . then qualitatively discuss other mod¬ 
els (E02)/ 


5.2.1. High Outflow Velocity 

An outflow velocity larger than V/xp ~ 300 km s -1 
can reduce AuLya.r because Lya photons would drop out 
of resonance with Hi atoms in the outflowing gas (e.g., 
IVerhamme et al.l 120061 l20l5l) . However, our results of 
Lya radiative transfer fitting in 94.3.21 show that all ob¬ 
jects have small V exp of 100 — 200 km s -1 . Combined 
with the findings in 94.5.1l that these 14 xp are consistent 
with the observed outflow velocities, Ar a b s , we conclude 
that the high outflow velocity hypothesis is unlikely. 


5.2.2. Low N m 


We examine the low Ahi hypothesis. Although 
it is difficult to directly measure iVni in LAEs from 
observations, we have inferred it using the expand¬ 
ing shell model 1 >14. 3. 2D . If we exclude the blue- 
bump objects from the sample, modeled Lya pro¬ 
files and parameters are all consistent with the ob¬ 
served Lya profiles and several fundamental observ¬ 
ables. Thus, we consider the derived neutral hydro¬ 
gen column density, Ahi, to be reliable as well. Fig¬ 
ure HD is a plot of A^Lya.r against log(Ani) for the non¬ 
blue bump objec t s. We add results from the lit e rature : 
Verhamme et al. (120081) . ISchaerer fe Verhammd (120081). 
Vanzella et al.l ( 201fll ). and iDessauges- Zavadsky et al.1 
(120101) . These authors have also utilized the model 

used in th is study for z ~ _ 3 LBGs with various 

EW(Lya) (IVerhamme et al.l 120081) . a strongly lensed 
LBG with Lya absorpt ion at z ~ 2.73 (MS 1512-cB58) 
(ISchaerer fe Verhammd 120081) . a peculiar z = 5.56 [Niv 


emitter with EW(Lya) = 89 A (I Vanzella et al] 1201 04. and 
a lensed LBG with Lya absorp tion, “the 8 o’clock arc” 
(IDessauges- Zavadsky et al l 120101). We also add th e re- 
sults of lKulas et al.l ( 20121) and lChonis et all (120131) . al¬ 
though models used in these studies are different from 
the one used in this study. In the figure, objects with 
EW(Lya) > 30 A are colored in red and labeled as LAEs, 
while those with EW(Lya) < 30 A are colored in blue 
and labeled as LBGs. 

We caution readers that all the data points and error 
bars in Figure [Tl] are obtained assuming uniform expand¬ 
ing shell models (see 1 14.21 for how we have obtained the 
error bars of our data points). The results can be sig¬ 
nificantly changed once if we consider other theoretical 
models such as clumpy or patchy models (see 95.2.311 . 

The figure shows a clear cor relatio n between log(lVHi) 
and AuLya, r - As described in 94.3.21 the mean log(Ani) 
in 2 ~ 2 LAEs is log(Ani) = 18.9 cm -2 , which is more 
than one order of magnitude lower than those of z > 3 
LBGs, ~ 20.0. We have excluded the 5 blue bump ob¬ 


jects in our 12 LAEs in Figure fill for a secure discussion. 
However, we note that they have comparable TVhi values 
to the non-blue bump objects, and are consistent with 
the correlation. We conclude that the small Ai>Ly Q ,r in 
the LAEs can be well explained by the low Ahi hypoth¬ 
esis. 


5.2.3. Other Models 

Throughout the paper, we have as sumed a uniform 
expan ding s hell models con s truct ed bv IVerhamme et all 
(j2006l) and ISchaerer et al.l (|2011l) . In this section, we 
qualitatively discuss alternative mo dels such as clumpy 
and patchy models. IHansen fc Ohl (j2006| ) have analyt¬ 
ically investigated the Lya radiative transfer in multi¬ 
phase media, especially, the one through dusty opti¬ 
cally thick gas clumps. They have shown that radiative 
transfer strongly depends on covering factor, f r (see also 
IDiikstra fe Kramerll2012| : lLaursen et al.l I2013D. As can 
be seen from Figure 20 of IHansen fc Ohl (120061 ). asym¬ 
metric Lya profiles with small Ai>Lya, r can be repro¬ 
duced by clumpy models with low / c , i.e., a small av¬ 
erage number of interaction for L ya photons before es ¬ 
capi ng from the galaxy. Re cently, IBehrens et al.l (|2014D 
and IVerhamme et al.l (|2015t) have investigated the Lya 
radiative transfer in a shell with holes and/or cavities, 
i.e., CF < 1. These studies have shown that, if a shell 
has holes, the modeled Lya line profile through a trans¬ 
parent line of sight has AuLya.r = 0 km s -1 even if 
convolved with spectral resolut ions used for the obser - 
vations (see Figures 4 and 5 in IVerhamme et aLll2015D . 
We qualitatively test the hypothesis by comparing Lya 
and nebular emission line profiles. In the case of a patchy 
ISM, if we observe the galaxy through a transparent line 
of sight, observed Lya line profiles would be indistin¬ 
guishable from nebular emission line profiles. This is 
because the main Lya component is dominant and is 
not affected by the radiative transfer effect. As can be 
seen in Figure 1, COSMOS-08357 and COSMOS-43982 
have indistinguishable Lya and nebular line profiles. 
We conclude that at least two objects, COSMOS-08357 
and COSMOS-43982, could be explained by patchy ISM 
models. 

Thus, we have demonstrated that the small AuLya.r 
in LAEs can be also well reproduced by clumpy and/or 
patchy ISM. However, future detailed Lya modeling as¬ 
suming clumpy and patchy ISM are needed for a more 
definitive conclusion. Whichever hypothesis is the most 
relevant one, the key for the small Ai>Ly a ,r in LAEs would 
be the reduced number of resonant scattering of Lya pho¬ 
tons. 

Hereafter in this section, we focus on the results ob¬ 
tained from uniform expanding shell models with low 

N m . 

5.3. Interpretation of Low TVhi in LAEs 

We have shown that, on the assumption of uniform 
shell models, the most likely situation for the smaller 
AuLya.r in the present sample is that LAEs have low 
-ZVhi. We can envisage three possible scenarios for LAEs 
having low Ahi- 

First, it is possible that LAEs have a low Hi gas mass. 
Indeed, iPardv et al.l (I2014D have detected an Hi 21cm 
line for 14 local galaxies with Lya emission (Lyman Al¬ 
pha Reference Sample; lOstlin et al.1120141) . to find that 
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Fig. 11 . log(A r H I) plotted against Ai;i, y „ r. Red (blue) symbols correspond to objects with EW(Lya) p i 10 to larger (smaller) than 30 A. 
Filled circles and diamonds are the non blue bump LAEs obtained by MagE and LIRS, respectively. Filled squares show z ~ 3 objects given 
bv IVerhamme et~a l. (2008). A filled triangle denotes a z = 5.56 [Niv emit ter with EW(Lyct) = 89A|Vaiizclla ct ah 2 0101 ). An open triangle 
is a lensed LBG wit h Lyo absorption, cB58 (Schaerer & Vcrhammc 2008), while an invert ed triangle is a lens ed L BG with Lyev absorp tion, 
“the 8 o’clock arc” iDc ssaugcs-Zavadskv ct al. 2010). In addition, the objects studied in IKulas et all 120121 ) and ICh onis ct al. (20 131) are 
plotted as a circle and three squares with a cross inside, respectively. For the purpose of display, three LAEs given bv IGhonis et al.l 12013 ) 
which have similar AuLya.r and log(A r ni) values of 175 km s' 1 and 18 cm' 2 , are shifted toward x-axis. We stress that all the data 
points and error bars are obtained assuming uniform expanding shells. The positions and/or error bars of data points in the figure can be 
significantly changed if we consider other models. See the text ( 45.2.211 for the detail. 


the derived Hi gas mass tentatively anti-correlates with 
EW(Lyct). This trend is also consistent with theoretical 
predictions (private communications with T. Garel and 
C. Lagos). 

Second, if a galaxy has a high gas ionization state, 
ionizing photons would efficiently ionize the neutral gas 
in the ISM. This would decrease the thickness of the 
Hi gas in the outflowing shell, and lower their Nm- 
This picture is consisten t with the recent finding of 
iNakaiima fc Ouchi ( 2014 1) that LAEs have a significantly 
higher ionization state than that of LBGs at the same 
redshift. The high ionization state in L AEs wo ul d be 
due to their young stell ar populations (e.g.. lPirzkal et all 
l2007t lOno et all 1 201 Oli fal) . Young O- and B type stars in 
LAEs would efficiently produce ionizing photons, and re¬ 


duce the Ah i of the surrounding gas. 

Finally, in the case of a face-on galaxy, we would see 
a lower Ahi because Lyct photons woul d experience a 
short e r path length out of the d isk (e.g.. lVerhamme et al ' 


l2012t iZheng fc Wallacel 120141) . Indeec Pshibuva et al 


(12014a!) have statistically examined the cllipticity, an in¬ 
dicator of the inclination, for z ~ 2 LAEs using HST 
data. A weak trend has been found that high EW(Lya) 
objects are less inclined. 

A combination of these effects would reduce the num¬ 
ber of resonant scatterings of Lya photons. This would, 
in turn, decrease the Lya velocity offset, Avi vn r (e.g., 
McLinden et al.l 1201 It lHashimoto et al.l 120131 Erb et ali 


20141) and the Lya sp atial offset, dJiang et al .11201.'ll 


Shibuva et al.ll2014al) . 
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5.4. Implication of Small AuL y «,r and Aupeak in LAEs 

Recent theoretical studies (IBehrens et al.1 12014 
lYerhamme et al.ll2015T ) have proposed that the Ly a line 
profile can be used as a probe of Lyman continuum 
(LyC; A < 912A) leaking galaxies (LyC leakers). LyC 
leakers are thought to have contributed to cosmic reion¬ 
ization. Observationally, detections of LyC emission 
are cl aimed for LAEs and LBGs both spectroscopically 
(e.g.. iShaplev et~aLl 120061) and photom etrically (e.g., 
Ilwata et al.l 120091 and lNestor et al.lI2013T) . However, the 
success rate is very low possi bly because LyC le akers are 
e xtremely faint objects (e .g. JOuchi et al.ll2008fl . 

IVerhamme et ~ni (f2oi5h have investigated two scenar¬ 
ios for the ionizing photon escape: (1) the density 
bounded Hll regions with an extremely low IVhi value, 
log(AHi) < 17.0 cm -2 , corresponding to an unity op¬ 
tical depth for ionizing photons, (2) or a galaxy has a 
partial spatial covering fraction of the gas. They have 
shown that AuL ya , r (Av pea k, if the blue bump exists) is 
extremely small in these cases, Ai>L ya ,r 1$ 100 km s _1 
(Afneak 5, 150 km s” 1 ) (see also iJaskot &; Oevl 12014 
iMartin et al.ll2015l) . 

On the other hand, as described in H3.41 objects with 
a large EW(Lya) value tend to have a small AuL ya ,r 
(Aup ea k), implying that they are good candidates of 
LyC leakers. If the ionizing photon escape fraction, 
/ osc , is high, the rec ombination Lya line mig ht be 
weaken. However, as iNakaiima fe Ouchil (12014D have 
shown, EW(Lya) can be as high as 100 A even at 
/esc ~ 0.5. 

Thus, we propose that selecting objects with EW(Lya) 
as large as 100A is a promising way to search for LyC 
leaking galaxies. The merit of this candidate selection 
technique is that it does not require spectroscopy. 

6 . SUMMARY AND CONCLUSION 

We have presented the results of a Lya profile anal¬ 
ysis of twelve LAEs at z ~ 2.2 for which high spec¬ 
tral re solut ion Lya lines are obtain ed in lHashimoto et al.l 
(|2013l) and lShibuva et all (|2014bD with Magellan/MagE 
or Keck/LRIS. Two objects have been observed with 
both spectrographs. All twelve objects have detections 
of nebula emission lines which are used not only to de¬ 
fine the systemic redshift but also to infer the intrinsic 
FWHM of the Lya line. We have also derived the galac¬ 
tic outflow velocity from blueshifted low-ionization state 
(LIS) metal absorption lines with respect to the systemic 
redshift for three individual LRIS spectra as well as for a 
stacked spectrum of four MagE spectra. In addition, we 
have obtained stellar dust extinction from SED fit. The 
high spectral resolution Lya data in conjunction with 
these measurements have enabled us to perform detailed 
comparisons between observed and modeled Lya lines. 
Our main results are as follows. 

• We find that all 12 objects have Lya profiles with 
a main peak redward of the systemic redshift and 
five objects (six spectra) have a weak, secondary 
peak blueward of the systemic redshift (the blue 
bump). For a sample of 17 objects from our study 
and the literature with a resolved Lya line, we 
estimate the ratio of LAEs with a blue bump to 
be ~ 50%, which is slightly higher than that of 


LBGs. We have obtained AuL ya , r = 174 ± 19 km 
s -1 (AwL y a,b = — 316±45 km s _1 ), which is smaller 
than (comparable to) that of LBGs, AuL ya ,r — 400 
km s _1 (AuL y a,b = —367 ± 46 km s -1 ). 

• The high spectral resolution and sensitivity of 
Subaru/FMOS have enabled us to detect two- 
component [Om] profiles in two LAEs for the first 
time. While its origin is not clear, we find that even 
the FWHM of the broad component is as small as 
70 — 80 km s -1 . This excludes the possibility of its 
origin being AGN activity or powerful hot outflows. 

• We have applied th e uniform expanding shell 
model constructed by IVerhamme et al.l (|2006t) and 
iSchaerer et al.l (|2011l) to our sample. The model 
successfully reproduces not only Lya profiles but 
also the galactic outflow velocity measured from 
LIS absorption lines and the FWHM of nebu¬ 
lar emission lines for the non blue-bump objects. 
However, for the blue-bump objects, the intrinsic 
FWHMs of Lya predicted by the model is signifi¬ 
cantly larger than the observed FWHMs of nebular 
emission lines. 

• For the blue bump objects, we have tried another fit 
fixing the intrinsic FWHM of Lya to the observed 
FWHM of nebular emission lines. The position and 
flux of the blue bump are poorly reproduced. 

• To understand the large discrepancy between 
FWHMj nt (Lya) and FWHM(neb) in the blue- 
bump objects, we have examined if objects with 
and without a blue bump have different properties 
such as the Lya luminosity, stellar mass, and the 
merger fraction. We find no significant difference 
between the two samples. We propose that tak¬ 
ing into account Lya photons produced by gravi¬ 
tational cooling might simultaneously explain the 
large FWHMj nt (Lya) and the existence of a blue 
bump. 

• We quantitatively demonstrate that the small 
AwL y o:,r in LAEs can be well explained by uniform 
expanding shell models with neutral hydrogen col¬ 
umn density as low as log (IVhi) = 18-9 cm~ 2 . This 
value is more than one order of magnitude lower 
than that of LBGs, and is consistent with the recent 
findings that LAEs have a high ionization param¬ 
eter and a low Hi gas mass. These results imply 
that low IVhi is the key for the small AwL ya ,r as 
well as the Lya escape mechanism. However, we 
caution readers that our results are based only on 
uniform expanding shell models, and that future 
detailed modeling with clumpy and/or patchy ISM 
are needed for a definitive conclusion. 

• As an implication of the small AuL ya ,r and low IVhi 
in high EW(Lya) objects, we propose that target¬ 
ing high EW(Lya) objects would be an efficient 
way to search for Lyman Continuum leaking galax¬ 
ies from photometry data alone. 
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APPENDIX 

Figures 12 - 14 show 2D y 2 contours for CDFS-3865. All the grids are shown in grey dots. The blue (red) grids in 
these figures show those satisfying Ay 2 < 6.17 (11.8) above the raw minimum y 2 designated by the white dots, i.e., 
the 3 (5) a uncertainty in the parameters dPress et al.lll992H . 
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Fig. 12.— Upper and lower five panels show 2D \ 2 contours for \4xp and log(A^Hi), respectively, for CDFS-3865. The grids colored with 
blue (red) denote those within the 3 (5) a level from the minimum x 2 grid shown as a white dot. 






























































Fig. 13.— Upper and lower five panels show 2D x 2 contours for r a and 6, respectively, for CDFS-3865. The meaning of the colors is the 
same as in Figure 12. 
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